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ABSTRACT 

Mergers of two carbon-oxygen (CO) WDs have been considered as progenitors of Type la supernovae (SNe 
la). Based on smoothed particle hydrodynamics (SPH) simulations, previous studies claimed that mergers 
of CO WDs lead to an SN la explosion either in the dynamical merger phase or stationary rotating merger 
remnant phase. However, the mass range of CO WDs that lead to an SN la has not been clearly identified 
yet. In the present work, we perform systematic SPH merger simulations for the WD masses ranging from 
0.5 Mq to 1.1 Mq with higher resolutions than the previous systematic surveys and examine whether or not 
carbon burning occurs dynamically or quiescently in each phase. We further study the possibility of SN la 
explosion and estimate the mass range of CO WDs that lead to an SN la. We found that when the both WDs are 
massive, i.e., in the mass range of 0.9 M 0 <Mi 2 < 1.1 Mq, they can explode as an SN la in the merger phase. 

On the other hand, when the more massive WD is in the range of 0.7 Mq<Mi<Q.9 Mq and the total mass 
exceeds 1.38 Mq, they can finally explode in the stationary rotating merger remnant phase. We estimate the 
contribution of CO WD mergers to the entire SN la rate in our galaxy to be of ^ 9%. So, it might be difficult 
to explain all galactic SNe la by CO WD mergers. 

Subject headings: binaries: close — galaxies: evolution — supernovae; general — white dwarfs — hydrody¬ 
namics 


1. INTRODUCTION 

SNe la play the important roles in the determination of 
cosmological parameters as luminous standard candles (e.g., 
Riess et al. 1998; Perlmutter et al. 1999) and in the chemi¬ 
cal evolution of galaxies as major sources of iron group el¬ 
ements (e.g., Kobayashi et al. 1998). However, their pro¬ 
genitors are still not identified yet (e.g., Maoz et al. 2014). 
They are considered as a thermonuclear explosion of a CO 
WDs in a binary system, in which the WD accretes mass 
from its companion and the WD mass approaches the Chan¬ 
drasekhar mass (Mch 1.4 Mq). But, it is still controversial 
whether its companion is a non-degenerate star, i.e., single 
degenerate (SD) model, or a degenerate star, i.e., double de¬ 
generate (DD) model. In the SD model, a CO WD accretes 
hydrogen/helium-rich gas from the companion and increases 
its mass upto Mca- Finally, carbon burning starts at the cen¬ 
ter of the CO WD and explodes as an SN la (Whelan & Iben 
1973; Nomoto 1982; Hachisu et al. 1996, 1999a,b). On the 
other hand, in the DD model, both the components are CO 
WDs. Because binaries lose their orbital angular momentum 
by emitting gravitational waves, they will eventually merge. 
If their total mass exceeds Mch, the binary finally explodes as 
an SN la (Iben & Tutukov 1984; Webbink 1984). 

Some clues on the progenitors of SNe la were found in re¬ 
cent observations. In particular, neither surviving compan¬ 
ions nor signatures of them were detected in some SNe la 
(e.g., Schaefer & Pagnotta 2012). This fact supports the DD 
system, although no detection of companions can also be ex- 
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plained by the SD model (Di stefano et al. 2011; Justhum 
2011; Hachisu et al. 2012). In some SNe la, signatures of 
circumstellar matter (CSM) were detected (e.g., Patat et al. 
2007; Blondin et al. 2009; Simon et al. 2009; Sternberg et al. 

2011). This supports the SD model (e.g., Maoz et al. 2014, 
for a recent review), although this can be explained by some 
DD models (Raskin & Kasen 2013; Shen et al. 2013; Soker 
et al. 2013). Observations do still not clarify which model is 
the main progenitors of SNe la, i.e., SD or DD (or other type) 
system. 

The DD model has a theoretical difficulty as a progenitor 
model. Some theoretical studies indicated that CO WD merg¬ 
ers can not become an SN la, but collapse to a neutron star 
(e.g., Saio & Nomoto 1985, 2004; Nomoto & Kondo 1991). 
These studies calculated only the evolution of merger rem¬ 
nants after DD systems merged. Their calculations were one¬ 
dimensional (ID) spherically symmetric and assumed station¬ 
ary state. However, the merger of a DD system is a three- 
dimensional (3D) and dynamical event, so multi-dimensional 
hydrodynamical simulations are necessary to reach a definite 
conclusion. Benz et al. (1990) used 3D smoothed particle hy¬ 
drodynamics (SPH) code and simulated mergers of two CO 
WDs. Following this work, there were several studies on 
mergers of DD systems (Rasio & Shapiro 1995; Segretain & 
Chabrier 1997; GueiTero et al. 2004; Yoon et al. 2007; Loren- 
Aguilar et al. 2009; Fryer et al. 2010; Pakmor et al. 2010; Dan 
et al. 2011; Raskin et al. 2012, 2014; Moll et al. 2014). These 
studies concluded that some DD pairs can explode as SNe 
la. Such successful models can be divided by the dynamical 
phase when the SN la explosion occurs. 
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Pakmor et al. (2010) simulated mergers of massive CO 
WDs (~0.9 Mq) and found that carbon detonation initiates 
during the dynamical merger phase because of the compres- 
sional heating by the disrupted secondary which violently ac¬ 
cretes onto the primary. The binary system finally explodes 
as a sub luminous SN la. Pakmor et al. (2012a) also simulated 
mergers of more massive CO WDs (1.1 -1-0.9 Mq) and found 
that the system leads to a normal SN la. 

If carbon detonation does not initiate in the dynamical 
merger phase, the remnant undergoes three phases. The first 
phase is the early remnant phase (e.g., Shen et al. 2012; 
Kashyap et al. 2015), 100-1000 s after the secondary is com¬ 
pletely disrupted. In this phase, the merger remnant does 
not reach a quasi-stationary state yet and still has small non- 
axisymmetric structures. The second phase is the viscous evo¬ 
lution phase (Schwab et al. 2012; Shen et al. 2012; Ji et al. 
2013), 10^-10* s after merging. In this phase, the remnant 
reaches a quasi-stationary, axisymmetric state and it evolves 
in a viscous timescale. The third is the thermal evolution 
phase (Saio & Nomoto 1985, 2004; Yoon et al. 2007; Shen 
et al. 2012), ^ 10^ years. If off-center carbon burning occurs 
in these three phases before the rotating core of the remnant 
reaches Mch, it likely converts the CO WD to an oxygen- 
neon-magnesium (ONeMg) WD. The ONeMg WD finally 
collapses to a neutron star when its core mass reaches Mca- 
Yoon et al. (2007) simulated the merger of CO WDs with 
masses of 0.9 -I- 0.6 Mq and followed the evolution of the 
merger remnant. They found that the remnant can avoid off- 
center carbon burning and explodes as an SN la in the thermal 
evolution phase, if it satisfies several conditions. 

Thus, these works concluded that some DD systems can be¬ 
come SNe la. However, the mass range of CO WDs that lead 
to SNe la is not clarified yet. In this work, we simulate merg¬ 
ers of CO WDs with a mass range of 0.5-1.1 Mq with our 
SPH code until the end of the early remnant phase, and iden¬ 
tify the mass range of CO WDs that lead to SNe la. Using the 
obtained mass range of CO WDs, we estimate the contribu¬ 
tion of CO WD mergers to the entire SNe la in our galaxy. 

Similar parameter surveys have already been done by Dan 
et al. (2012, 2014) and Zhu et al. (2013). Zhu et al. (2013) 
mainly focused on the (early) remnant phase and the resolu¬ 
tions of their simulations are lower than our study. Although 
Dan et al. (2012,2014) covered the merger phase and remnant 
phase, their numerical resolution is much lower than ours. We 
expect that SN la explosions occur not only in the early rem¬ 
nant phase but also in the dynamical merger phase. In the 
present work, we adopt higher resolution than two times those 
of the previous works (e.g., Zhu et al. 2013). This is because 
the numerical resolution is one of the most important param¬ 
eters to identify the initiation of detonation in the dynamical 
mergerphase (Pakmor et al. 2012b). Therefore, we adopt four 
different resolutions and check the numerical convergence. 

To investigate the possibility that carbon burning leads to an 
SN la, we check the density and temperature of SPH particles 
and identify carbon burning by the condition of tqc < "TcooI, 
where tcc is the timescale of carbon burning and TcooI is 
the cooling timescale. In the dynamical merger phase, TcooI 
is the dynamical timescale of adiabatic expansion. On the 
other hand, in the (early) remnant phase, TcooI is the timescale 
of neutrino cooling. These are necessary but not sufficient 
conditions for a thermonuclear explosion. In this sense, our 
results would not be conclusive but are sufficiently sugges¬ 
tive. Since the temperature has numerical noise in our SPH 


simulation, we have to treat the temperature carefully. This 
noise comes from fluctuations of density and internal energy, 
which arise due to numerical resolution, finite neighbor par¬ 
ticles, and the accuracy of time-integration. We use both the 
raw and smoothed temperatures to estimate the noise in our 
SPH simulation (Dan et al. 2012). The raw temperature is 
the temperature that our SPH code originally generates. The 
smoothed temperature is an averaged temperature calculated 
from neighbor particles, which could avoid large numerical 
noises (defined in Section 3.1). In the present version of our 
SPH simulation, we do not include nuclear reactions because 
they are so sensitive to temperature noises and possibly en¬ 
hance them erroneously. 

This paper is organized as follows. Section 2 briefly de¬ 
scribes our numerical methods. In section 3, we show our 
results and then estimate the rate of SNe la coming from the 
DD merger systems and their contribution to the entire SNe 
la in our galaxy. In section 4, we compare our results with 
previous works and discuss the dependence of our numeri¬ 
cal results on the resolution. Finally we conclude the present 
work in Section 5. 

2. METHODS 

Here, we present a brief summary of our numerical method. 
The details of the numerical method are already described in 
Nakasato et al. (2012) and Tanikawa et al. (2015). 

2 .1. Numerical Code 

SPH is the Lagrangian mesh-free particle method devel¬ 
oped for the study of astrophysical fluid phenomena (Lucy 
et al. 1977; Gingold & Monaghan 1982). Nowadays, it is 
applied to more various subjects, e.g., engineering, meteorol¬ 
ogy, and Hollywood movies. Recent good reviews about SPH 
codes are available, e.g., in Monaghan (2005) and Rosswog 
(2009). Our formulation is the one called "vanilla ice" SPH 
formulation in Rosswog (2009). Our basic SPH equations 
consist of the equation of continuity, equation of motion, and 
energy equation for self-gravitating fluid in Lagrangian for¬ 
mulation. 

We use "OcTree On OpenCL"(OTOO) code for our 3D 
simulations of CO WD mergers, which is developed for var¬ 
ious particle simulations of astrophysical fluid phenomena 
(Nakasato et al. 2012). This code implements the octree 
method (Barnes & Hut 1986) to calculate gravity and neigh¬ 
bor particles for SPH. It is optimized for multiple CPUs and 
GPUs in heterogeneous computational resources. 

In this code, we adopt the 3rd-order spline kernel and its 
derivative is modified in the same way as that proposed by 
Thomas & Couchman (1992) to avoid the pairing instability. 
The smoothing length is determined to keep the average num¬ 
ber of neighbors being about 75 in every step (Thacker et al. 
2000). The formulation of the artificial viscosity is basically 
the same as Monaghan (1992), but the viscosity coefficient 
is time-dependent (Morris & Monaghan 1997). More detail 
explanation can be seen in Rosswog et al. (2000) and Mon¬ 
aghan (2005). We also introduce the Balsara switch to shut 
off the artificial viscosity in no shock regions (Balsara 1995). 
We adopt a leap-frog scheme for time integration. 

We use the Helmholtz equation of state (EOS) (Timmes & 
Swesty 2000) and assume a uniform chemical composition of 
50% carbon and 50% oxygen. This chemical composition is 
not changed throughout the merger simulation because we do 
not include nuclear reactions. 
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2 .2. Initial Setup 

The initial setup of our simulations is similar to those of 
Rasio & Shapiro (1995) and Dan et al. (2011). We separately 
generate each CO WD from spherically symmetric density 
profiles of a perfectly degenerate star, and set the tempera¬ 
ture at 10^ K everywhere. To reduce numerical noise in the 
mapping from the ID spherically symmetric density profile to 
our 3D SPH density distribution, we relax the SPH particles 
for 20 physical seconds with velocity damping force but with¬ 
out the evolution of internal energy. For the damping force, 
we adopt 


dt Cdamp^^ 

where v, is the velocity of /th particle, Cdamp is the inverse of 
relaxation timescale and we fix Cdamp = 128.0. 

Next, we put two CO WDs to the distance enough to avoid 
the Roche lobe over flow (RLOF). We assume that they are on 
a circular orbit and their spins synchronize each other with the 
orbital motion. To finish setting the initial condition, we grad¬ 
ually decrease the separation and relax them with the above 
damping force and the evolution of internal energy in a co¬ 
rotating frame. When an SPH (;th) particle of the secondary 
approaches closely enough the LI point, i.e., 

Ilrsecd-rull <0.2R2, (2) 

we stop the calculation. Here Tsec,; is the position of /th par¬ 
ticle of the secondary, Tli is that of LI point, and R 2 is the 
effective radius of the secondary. 

Finally, we transfer the SPH particles from the co-rotating 
frame to a rest frame. 

3. RESULTS 

We summarize the results of our simulations in Table AL 
Because we aim to identify the mass range of CO WDs that 
lead to an SN la, our initial models cover the entire mass range 
of CO WDs, i.e., fromMwD = 0.5 to 1.1 Mq. We prepare bi¬ 
nary models of 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, and 1.1 Mq, and 
perform the simulations of mergers for all 28 mass combina¬ 
tions. Since we also investigate the dependence of our results 
on the numerical resolution, we perform the same simulations 
with different resolutions, which are lOA:, 50A:, lOOA:, 500A: 
SPH particles per one solar mass (here A: = 1,024). All sim¬ 
ulations were performed from the start of RLOF until the for¬ 
mation of almost stationary remnant, several hundred seconds 
after the secondary is completely disrupted. We divide the dy¬ 
namical evolution of a merger into the following two phases. 
(1) The merger phase is the period during which the secondary 
dynamically accretes onto the primary. (2) The (early) rem¬ 
nant phase is from the complete disruption of the secondary 
to the stage that the system reaches almost a quasi-stationary 
state. 

Figure 1 shows the density profiles in the equatorial plane 
for an example of our simulations, whose mass combination 
is 1.1 -t 0.9 Mq and resolution is 500A: Mq . The merger phase 
covers the first (t = 25 s) to fifth (f = 135 s) panels, and 
the early remnant phase corresponds to the sixth (f = 240 s) 
panel. Their morphological structures are consistent with the 
previous works (e.g., Pakmor et al. 2012a). 

First we check whether dynamical carbon burning starts or 
not in the merger phase for all simulations because it is a nec¬ 
essary condition for an SN la explosion. Then, we check 
steady carbon burning in the early remnant phase. If it oc¬ 
curs in the early remnant phase, carbon burning converts a 


CO WD into an ONeMg WD and the merger remnant will not 
become an SN la. Although we perform simulations with four 
different resolutions, we first focus on the highest resolutions 
(= 500k Mg'). 


3.1. Merger Phase 

Pakmor et al. (2010) first suggested that CO WD mergers 
lead to an SN la in the merger phase. They called such a 
model the (carbon-ignited) violent merger (VM) scenario. In 
this scenario, matter of the secondary violently accretes onto 
the primary and such violent accretion causes dynamical car¬ 
bon burning. As a result, detonation wave would be formed 
and propagate into the primary with converting its carbon- 
oxygen into iron group elements. Finally, the system explodes 
as an SN la. Because our simulations cannot directly resolve 
the initiation of detonation, we try to judge the occurrence of 
dynamical carbon burning in the merger phase. This condition 
is a necessary condition for an SN la explosion. 

For this purpose, we extract the highest temperature particle 
in the merger phase for all simulations. The condition for 
dynamical carbon burning is 


TCC < Tdyn, 

where tcc is a carbon burning timescale defined by 

CpT 

Tcc = -) 

ecc 


(3) 

(4) 


and Tdyn is a dynamical timescale (Nomoto 1982) defined by 


where Cp is the specific heat at constant pressure, ecc is the 
energy generation rate of carbon burning. We calculate the 
both timescales for each particle in the merger phase and ex¬ 
amine whether the particles satisfy Equation (3). In this work, 
Cp is derived from the Helmholtz EOS of Timmes & Swesty 
(2000) and ecc is the same as that of Dan et al. (2014), origi¬ 
nally proposed by Blinnikov (& Khokhlov (1987). Its formu¬ 
lation is 

ecc = pqcAT Yc exp)-^/?’^/^ +/cc), (6) 


where qc = 4.48 x 10'* erg mol * (Blinnikov & Khokhlov 
1987), At = 8.54xl0267’9^/V”^/^s-'mor'cm*, 

Ti) = r/10®K, Tga = Tg/(l+0.067Tg), Q = 84.165 (Eowler 
et al. 1975). The carbon abundance is calculated as 
Yc = «c/(piVa) = 0.033 mol g“*, where nc is the number den¬ 
sity of carbon and N^. is the Avogadro number. A screening 
factor fee is ignored here, because we focus on the start of 
dynamical carbon burning and the factor of self-acceleration 
for nuclear burning described in Erank-Kamenetskii (1967) is 
not applied to the initiation of carbon burning. 

Eigure 2 shows the density and temperature of a particle 
with the highest temperature in the merger phase for all mass 
combinations (their resolutions are 500k Mg), similar to Eig¬ 
ure 12 of Dan et al. (2014). The shapes and colors of sym¬ 
bols indicate the primary’s mass and the total mass, respec¬ 
tively. Solid lines indicate rcc = fiiyn, and dashed lines do 
Tcc = 0. iTdyn. Eor the mass combinations above the solid line, 
dynamical carbon burning occurs in the merger phase. So, the 
merger of CO WDs would lead to an SN la explosion. 

In our SPH simulation, physical raw temperature of each 
particle has numerical noise. So we adopt another definition 
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Fig. 1.— Density profiles in the equatorial plane for the dynamical evolution 
resolution is 500k Mq . Colors indicate density in a logarithmic scale. 

of temperature to reduce the effect of noise, i.e., smoothed 
temperature. It is defined by 

Z,i = J2—TjW(rij,h,j), (7) 

where mj, pj, Tj are the mass, density, temperature of jth par¬ 
ticle, respectively. = ||r, -rj||, and hjj = {hi + hj)/2 is the 
average of the smoothing lengths of ith and y'th particles. 

Figure 2(a) and 2(b) show the results for raw temperature 
and smoothed temperature thus defined, respectively. It is 
clear that all symbols in Figure 2(a) move to a lower place 
in Figure 2(b). As a result, the number of mass combina¬ 
tions above the solid line of rcc = rdyn decreases. Figure 2(b) 
of smoothed temperature shows that the mass combinations 
above the solid line would certainly trigger dynamical carbon 
burning in the merger phase. 

3.2. Remnant Phase 

If dynamical carbon burning does not occur in the merger 
phase, the merged object goes into a stationary phase, i.e., 
the (early) remnant phase. Figure 3(a) and 3(b) show the 
density and temperature, respectively, of the merger remnant 
whose mass combination and resolution are 1.1 -1-0.9 Mq and 
500k Mq. The remnant consists of three components, i.e., a 
cold core, hot envelope, and outer disk. The structure of such 
a remnant has been studied in the several works as already 
mentioned in Section 1. Although there are small differences 
among these previous works, their results are almost consis¬ 
tent with each other. Our results are also consistent with these 
previous results. 

It has long been discussed that such a hot envelope gradu¬ 
ally accrete onto a cold core because of angular momentum 
loss by some mechanisms (e.g. viscosity or magnetic field). 

If off-center carbon burning occurs quiescently during such 
an accretion phase, carbon deflagration waves propagate into 
the core. At last, the whole core is converted into an ONeMg 
WD (e.g., Saio & Nomoto 1985, 1998, 2004). In such a case. 


of our merger simulation. The mass combination is 1.1 +0.9 Mq, and the 

it can not explode as an SN la, even if its total mass exceeds 
Men- Instead, it collapses to a neutron star (Nomoto & Kondo 
1991). On the other hand, if off-center carbon burning does 
not occur, the core remains unchanged as a CO WD and sur¬ 
rounding matter continues to accrete onto the core. When 
the core mass exceeds Men, it explodes as an SN la. It is 
critically important to examine whether carbon burning starts 
quiescently off-center in the remnant phase. 

We examine, in the same way as that of dynamical carbon 
burning in the merger phase, whether or not off-center carbon 
burning occurs in the remnant phase. Carbon burning quies¬ 
cently occurs near the boundary between the cold core and hot 
envelope, if the condition for carbon burning, i.e., 

fee < fi/5 (8) 

is satisfied in the remnant phase. Flere, is a timescale 
of neutrino cooling and we use the description of Itoh et al. 
(1996) for calculating it from the density and temperature of 
SPH particles. We find the highest temperature particle in 
the remnant phase for all mass combinations, and examine 
Equation (8). If the highest temperature particle satisfies the 
condition, we regard that off-center carbon burning starts and 
converts the CO core into an ONeMg core. Then, the system 
finally collapses to a neutron star if the total mass exceeds 
Mch. On the other hand, if there are no particles that satisfy 
the condition and the total mass of the remnant exceeds Mch, 
we consider that the remnant becomes an SN la. 

Figure 4 is the same plot as Figure 2, but for the remnant 
phase. Magenta solid lines indicate rcc = fi/- Figure 4(a) 
and 4(b) show results of the raw and smoothed temperatures, 
respectively. Off-center carbon burning occurs in the remnant 
phase for the mass combination models above the magenta 
solid line. They would finally collapse to a neutron star. On 
the other hand, the models below the line would become an 
SN la if the total mass exceeds the Chandrasekhar mass, i.e.. 
Ml + M 2 > 4/ch' 

We have to follow the viscous and thermal evolution phases 
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Fig. 2.— Density and temperature of the highest temperature particle in 
the merger phase for all mass combinations. Numerical resolutions of these 
models ai‘e 500k Mq . Colors of symbols indicate the total mass of the sys¬ 
tem as indicated in the rightside of the figures, and shapes of symbols do the 
mass of the primary. Filled squares are the 1.1 Mq primary, filled circles 
1.0 Mq, filled triangles 0.9 Mq, filled inverted triangles 0.8 Mq, filled di¬ 
amonds 0.7 Mq, open squares 0.6 Mq, open circles 0.5 Mq. Solid lines 
indicate tqc = '^dyn^ dashed lines do tqc = O-lr^yn. (a) Raw temperature 
of SPH particles, (b) Smoothed temperature, Ts, of SPH particles defined by 
Equation (7). 


for obtaining definite conclusion on the off-center carbon 
burning. However, we stop our SPH calculation at the end 
of the early remnant phase because our SPH code includes 
no physical viscosities (e.g., magnetic viscosity). In the vis¬ 
cous and thermal evolution phases, the hot envelope further 
accretes onto the cold core and compress itself on the cold 
core. The temperature and density might further increase. As 
a result, carbon ignites off-center even for the cases of no off- 
center burning in the early remnant phase. In fact, Shen et 
al. (2012) and Schwab et al. (2012) followed the evolution of 
merger remnants and showed that off-center carbon ignition 
starts in the viscous and thermal evolution phases for some 
cases. Yoon et al. (2007) performed SPH simulation of CO 
WD merger whose mass combination is 0.9 -I- 0.6 Mq and 
further followed the evolution of the merger remnant with a 
ID stellar evolution code. They found that off-center car¬ 
bon burning can be avoided when the highest temperature 
is lower than the threshold for carbon ignition, angular mo¬ 
mentum loss occurs with a timescale longer than that of neu¬ 
trino cooling, and the mass accretion rate is M < 5x10“^ 
to 10“^ Mq yr“'. Our present condition for off-center carbon 
burning is posed only for the (early) remnant phase but not ap¬ 
plied yet for the viscous and thermal accretion phases. In this 
sense, our results only for the early remnant phase are not def- 



Fig. 3.— Structure of a merger remnant at about 200 s (~7.6 Porb.init) af¬ 
ter the secondary is disrupted completely. Here, Porb.init is the initial orbital 
period. Its mass combination is 1.1 + 0.9 Mq and the numerical resolution is 
500A: Mg , which is the same model as in Figure I. (a) Density profile (in a 
logarithmic scale of g cm“'^) in the —z plane, (b) Temperature structure (in 
a linear scale of 10® K) in the x—z plane. The central part becomes slightly 
hot. This is caused by numerical noise and artificial viscosity, so it has no 
physical meaning. 

inite answers. We leave such viscous and thermal evolutions 
of the merger remnants in our future works. 

3.3. White Dwarf Mass Combinations For SNe la 

Now, we have obtained the mass range of CO WDs which 
possibly lead to an SN la. Using this mass range of CO 
WDs, we estimate their contribution to the entire SNe la in 
our galaxy. We consider four paths that CO WD mergers 
would follow. The first one is the VM path, the condition of 
which is that dynamical carbon burning occurs in the merger 
phase. The systems satisfying this condition would explode 
as an SN la immediately after merging (Pakmor et al. 2010, 
2011, 2012a). When dynamical carbon burning does not oc¬ 
cur in the merger phase, the system enters the remnant phase 
and the disrupted secondary surrounds the primary. If its total 
mass exceeds Mch and off-center carbon burning does not oc¬ 
cur during the accretion phase, carbon burning occurs at the 
center of the CO WD and it would finally explode as an SN 
la. We regard this evolutionary path as the accretion induced 
explosion (AIE) path. On the other hand, if off-center car¬ 
bon burning occurs, the core of the remnant will be converted 
into an ONeMg WD and then collapse to a neutron star when 
the core mass exceeds Mch. This is the accretion induced col¬ 
lapse (AIC) path (Saio & Nomoto 1985,1998, 2004; Nomoto 
& Kondo 1991). When dynamical carbon burning does not 
occur in the merger phase and the total mass of the system 
does not exceed Mch, the system would form a massive CO 
WD. We call this evolutionary path as the massive white dwarf 
(MWD) path. 

Figure 5 shows all mass combinations of our simulations 
with 500k Mq particles and identifies which path they take. 
Colors of symbols indicate the four paths, i.e., the VM (red), 
AIE (green), AIC (blue), and MWD (magenta) paths. Among 
these four paths, the VM and AIE paths are possible paths to 
SNe la. 
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Fig. 4.— Same as Figure 2, but for the remnant phase. Magenta solid lines 
indicate tqq = . 

3.4. SN la Rate 

We estimate the rate of SNe la which originate from CO 
WD mergers and their contribution to the entire SNe la in our 
galaxy. Here, we assume that all mergers of CO WDs sat¬ 
isfying the VM or AIE condition can explode as an SN la. 
Strictly speaking, for the VM path, dynamical carbon burning 
is a necessary, but not a sufficient condition for carbon deto¬ 
nation. So it is not trivial if those mergers explode as an SN la. 
On the other hand, we have to follow the evolution of merger 
remnants for a much longer time (e.g., Yoon et al. 2007) in 
order to identify their final fates, the AIE (exploding as an SN 
la) or AIC (collapsing to a neutron star) path. In this sense, 
our estimate is just an upper limit. 

The SN la rate depends on the Hubble type and stellar 
mass of a galaxy. According to Li et al. (2011), in SBc 
type galaxies with similar stellar mass to our galaxy, the SN 
la rate is about 1.1x10“'^ yr“* Mg. Badenes et al. (2012) 
estimated the merger rate of binary WDs in our galaxy as 
1.4x 10“'^ yr“' Mg . The mass distribution of binary WDs in 
our galaxy is still uncertain because there are small samples 
even in our neighborhood. We assume that both the primary 
and the secondary follow the mass distribution of single DA 
WDs in our galaxy derived from SDSS-DR7 (see, e.g.. Figure 
10 of Kleinman et al. 2013). Then we can calculate the rate of 
CO WD mergers that satisfy the condition of each scenario. 

The merger rate is about 1.4x10“'^ yr“'Mg 
(0.14 x 10-'5 yf'Mg') for the VM, 5.4x10-'^ yr'Mg' 
(8.9x10“'^ yr'^Mg) for the AIE path, and 
6.8x10“'^ yf'Mg' (9.0x10“'^ yV'Mg') for the both 


Fig. 5.— The outcome of our merger simulations for the (a) raw and (b) 
smoothed temperatures. Red symbols mean the VM path, green symbols 
the AIE path, blue the AIC, and magenta the MWD path. Two black lines 
indicate that the primary and secondary have the same mass and that the total 
mass equals M\ +M 2 = Mch '^1-4 Mq . Mass combinations of red and green 
symbols result in an SN la. 

paths, if we adopt the case of raw temperature (smoothed 
temperature). This is only ^ 9% of the entire galactic SNe la. 
Therefore, at least in our galaxy, DD merger systems might 
not dominate the progenitors of SNe la. 

4. DISCUSSION 

4.1. Comparison With Previous Studies 

We compare the present results with previous studies. In 
particular, we mainly focus on two mass combinations. One 
is 1.1 4-0.9 Mg and the other is 0.94-0.6 Mg, because these 
two combinations were well studied in the previous works. 

4.1.1. I.I-hO.9Mq 

First we compare our result of 1.14- 0.9 Mg and 100k Mg 
(i.e. the total number of SPH particles is about 2x 10^) with 
that of Pakmor et al. (2012b). Figure 6(a) shows the time 
evolutions of the orbital separation and Figure 6(b) shows the 
number of particles having temperature higher than 2 x 10“^ K. 
Our Figure 6 should be compared with Figure 4 and Table 1 
of Pakmor et al. (2012b). 

It should be noted that the WDs in our model merge more 
quickly (about 100 s) than those in Pakmor et al. (2012b) 
(about 600 s) although we start our simulation with the ini¬ 
tial condition similar to theirs. This is because our initial sep¬ 
aration (^1.67x10® cm) is less than Pakmor et al. (2012b) 
(^1.93x10^ cm). Indeed, Figure 6(a) and 6(b) resemble the 
case of their smaller initial separation. We suppose that the 
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Fig. 6.— Time evolution of (a) the orbital separation and (b) the number of 
high temperature (>2.0x 10^ K) particles. The mass combination of WDs is 
1.1 + 0.9 Mq and the numerical resolution is lOOt: Mg . 

main reason of the difference is the relaxation method of a 
single WD. As mentioned in Section 2.2, we relax a single 
WD with a velocity-dependent damping force but without the 
evolution of internal energy. On the other hand, Pakmor et al. 
(2012b) changed the damping timescale (see their Equation 
14). After the relaxation method of Pakmor et al. (2012b) is 
applied to our models, the radii of the relaxed WDs are a few 
percent larger than our original ones. Therefore, their separa¬ 
tions when the RLOF starts are larger than ours and it results 
in a longer merging time. 

If the initial separation is smaller when the RLOF starts, the 
mass transfer tends to occur more violently and the secondary 
is completely disrupted in a few orbital periods (Dan et al. 
2011). As a result, accreted matter is more strongly shock- 
heated and dynamical carbon burning easily occurs. Although 
Pakmor et al. (2012b) concluded that the initial condition is 
not so important for dynamical carbon burning, we should 
note that our highest temperatures in the merger phase could 
be overestimated. 

4.1.2. O.9 + O.6M0 

We also compare our result of 0.9 -I- 0.6 Mq and 100k Mq 
with the previous studies of Yoon et al. (2007), Dan et al. 
(2011), and Zhu et al. (2013). Since these studies have simi¬ 
lar resolution (a few x 10^ SPH particles in all), it is suitable 
for comparison. Figure 7 shows the time evolution of highest 
temperature in each density zone. Since no dynamical car¬ 
bon burning occurs (see Figure 5), we examine whether or not 
quiescent carbon burning occurs in the remnant phase, i.e., at 
t ~ 300 s in Figure 7. The temperature is about 6x 10* K and 
this is consistent with the above three previous studies. 

We compare our Figure 7 with Figure 4 in Yoon et al. 
(2007). For the highest temperature in the merger phase (at 
t ^ 100 s in Figure 7), their results of 1.7 x 10® K is higher than 
ours (1.3xl0®K), although dynamical carbon burning does 
not occur in the merger phase for both ours and theirs. We 
suppose that the difference in the highest temperature comes 
from the difference in the initial condition. Dan et al. (2011) 


N = 1 OOkMsun ) 0.9+0.6Msun 



Fig. 7.— Evolution of highest temperature in each density range, 
logp (g < 5.0 (red), 5.0 < logp < 5.5 (green), 5.5 < logp < 6.0 (blue), 
6.0 < logp < 6.5 (black), for the mass combination of 0.9+0.6 Mq and the 
resolution of 100<: Mg . 

reported that the morphology of merger remnant could be af¬ 
fected by the initial condition. They found that the remnant 
whose initial separation is larger has longer trailing arm than 
the one whose initial separation is smaller (see Figure 10 of 
Dan et al. 2011). On the other hand, we find that the high¬ 
est temperature in the remnant phase is barely affected by the 
initial condition (see Figure 7). Similar discussion appears in 
Tanikawa et al. (2015). 

Other factors might affect the results of our simulations. For 
example, Zhu et al. (2013) and Dan et al. (2014) performed 
merger simulations of non-spinning WDs and found that the 
structures of such merger remnants are different from those in 
mergers of synchronously spinning WDs, especially for the 
case of nearly equal masses. In non-spinning cases, the high 
temperature region is formed near the center of the merger 
remnant, and carbon burning might occur in that region. Since 
it is still uncertain whether binary WDs maintain synchroniza¬ 
tion until their merger, our results for the AIE path might be 
changed. The topic of synchronization is out of the scope in 
this paper, so we leave this in our future works. 

4.2. Numerical Resolution 

Pakmor et al. (2012b) concluded that the numerical reso¬ 
lution of simulation is one of the most important factors for 
carbon burning in the merger phase because very high reso¬ 
lution is required to identify very small hot spots. Therefore, 
we examine the dependence of our results on the numerical 
resolution. We perform the same simulations with four differ¬ 
ent resolutions, i.e., 10k, 50k, 100k, and 500k per one solar 
mass. The highest temperature is critical both for dynamical 
carbon burning and quiescent carbon burning, so we focus on 
the dependence of the highest temperature on the resolution. 

Figure 8(a) and 8(b) show the dependence of the highest 
temperature in the merger phase on the numerical resolution 
for the (a) raw and (b) smoothed temperatures, respectively. 
The highest temperature increases with the number of SPH 
particles. In other words, our simulation does not converge 
yet and the final fates of some models could be changed from 
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Fig. 8.— Dependence of highest temperature on the numerical resolution 
in the merger phase for the (a) raw and (b) smoothed temperatures. The 
horizontal axis is the number of SPH particles per a solar mass. The vertical 
axis is the highest temperature. Shapes and colors of symbols have the same 
meaning as those in Figure 2. The highest temperature tends to increase with 
the numerical resolution. 

the AIC path to the VM path. Such a tendency was also re¬ 
ported in Pakmor et al. (2012b). In this sense, we must further 
increase the numerical resolution to definitely identify the fate 
of merger products, at least, more than 500k Mg . Figure 9(a) 
and 9(b) show the dependence of the minimum rcc/Tayn ra¬ 
tio in the merger phase on the numerical resolution for the 
(a) raw and (b) smoothed temperatures, respectively. The ma¬ 
genta dashed line indicates Tcc/rdyn = 1 -0. It tends to decrease 
as the numerical resolution increases. This trend is consistent 
with that of the highest temperature. For smoothed temper¬ 
ature, only very massive pairs (both masses ^1.0 Mg) can 
ignite carbon dynamically. 

On the other hand. Figure 10(a) and 10(b) show the high¬ 
est temperature in the remnant phase for the (a) raw and (b) 
smoothed temperatures, respectively. Comparing with the re¬ 
sults in the merger phase (Figure 8), the highest temperature 
depends barely on the resolution. Especially, for the smoothed 
temperature in Figure 10(b), it converges for almost all the 
mass combinations. This tendency of weak dependence on 
the resolution was already reported in the previous studies 
(Raskin et al. 2012; Dan et al. 2014). 

5. CONCLUSIONS 

We have performed SPH simulations of CO WD merg¬ 
ers for the mass combinations of 0.5-1.1 Mg from the start 
of the RLOF to the formation of a quasi-stationary merger 
remnant, and examined whether carbon burning occurs ei¬ 
ther in the merger phase or remnant phase. Using the results 


Fig. 9.— Dependence of minimum rcc/'Tdyn ratio on the numerical resolu¬ 
tion in the merger phase for the (a) raw and (b) smoothed temperatures. The 
horizontal magenta dashed lines indicate tqc /"rdyn = 1.0. 

of SPH simulations, we have investigated the mass range of 
CO WDs that possibly lead to an SN la in the merger phase 
or remnant phase. We have obtained the mass range as fol¬ 
lows. When the primary and secondary are as massive as 
0.9 Mg<Mi_2<l.l Mg, the binary results in an SN la in 
the merger phase. On the other hand, when the primary is 
0.7 Mg < Ml < 0.9 Mg and the total mass of the binary ex¬ 
ceeds 1.4 Mg, they lead to an SN la in the remnant phase. 

From the obtained mass range, we have estimated the rate 
of SNe la coming from CO WD mergers in our galaxy. It is 
6.8 X10“'^ yr“' Mg if we use the results of our raw tempera¬ 
ture calculations, while it is 9.Ox 10“'^ yr“' Mg if we use that 
of the smoothed temperature. These are only less than 9% of 
the entire SN la rate. Therefore, it is unlikely that the mergers 
of CO WDs are the main progenitors of SNe la. 

Of course, the above estimate is not conclusive because 
of several uncertainties in our calculation. We have checked 
the dependence of the highest temperature on the numerical 
resolution both in the merger and remnant phases in order to 
examine the numerical convergence of our simulations. We 
have found that the highest temperature in the merger phase 
depends on the numerical resolution. It tends to increase with 
the resolution as already reported in Pakmor et al. (2012b). 
On the other hand, in the remnant phase, the highest tempera¬ 
ture depends barely on the numerical resolution. Therefore, it 
is necessary to increase the number of SPH particles, at least, 
up to >500k Mg , for definite conclusion. Additionally, our 
calculations for SNe la in the remnant phase is not sufficient. 
In order to obtain the decisive conclusion, we have to follow 
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the further evolution of the merger remnant, like Yoon et al. 
(2007). This is one of our future works. 
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Fig. 10.— Same as Figure 8, but for the highest temperature in the remnant 
phase. The highest temperature seems to converge. 


APPENDIX 

SUMMARY OF CO WD MERGER SIMULATIONS 

We summarize our numerical results in Table 1 and 2. The numeric values presented in the table are the primary mass. Mi, 
secondary mass, M 2 , initial orbital separation, ainit, raw and smoothed maximum temperatures, Tmax and Ts max. densities at 
the highest raw temperature and at the highest smoothed temperature, p{T„iax) and p{Tfi „iax), in the merger phase, Tmax.rem and 
7s,max,rem. p(7max,rem) and /^(Ts jjiax,rem) lit the remnant phase. 
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TABLE 1 

Summary of all calculated models 


Ml 

M2 


Tmax 

p(7max) 

^s,max 

p(7s,max) 

7max,rem 

p(7max,rem) 

^s,max,rem 

p(^s,max,rem) 

( Mq ) 

( Mq ) 

(10^ cm) 

(10^ K) 

(10*^ g cm-3) 

(10“’ K) 

(10'’ g cm-3) 

(10'* K) 

(lO*^ gcm-3) 

(lO** K) 

(lO*^ gcm-3) 







Resolution 

= lOA:M0-‘ 
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1.28 

0.28 

1.06 

0.5 

0.5 

2.35 

0.57 

1.06 

0.29 

0.65 

0.30 

0.76 

0.24 

0.49 






Resolution 

= 50A: M0-‘ 





1.1 

1.1 

1.23 

3.58 

6.13 

1.81 

8.88 

1.56 

3.02 

1.44 

3.09 

1.1 

1.0 

1.44 

3.10 

4.15 

1.67 

3.59 

1.27 

2.71 

1.14 

2.83 

1.1 

0.9 

1.65 

2.20 

2.87 

1.32 

2.66 

1.13 

2.94 

0.99 

2.36 

1.1 

0.8 

1.89 

1.94 

1.47 

1.16 

1.10 

1.08 

2.42 

0.93 

2.80 

1.1 

0.7 

2.13 

1.58 

0.39 

1.06 

3.90 

1.00 

1.96 

0.86 

1.76 

1.1 

0.6 

2.42 

1.52 

0.33 

0.88 

0.51 

0.94 

1.89 

0.82 

2.19 

1.1 

0.5 

2.79 

1.40 

0.67 

0.86 

0.60 

0.96 

1.40 

0.81 

1.49 

1.0 

1.0 

1.42 

2.77 

3.55 

1.31 

3.10 

1.11 

2.98 

0.99 

3.08 

1.0 

0.9 

1.62 

2.16 

2.67 

1.06 

3.01 

0.94 

1.87 

0.81 

2.73 

1.0 

0.8 

1.85 

1.93 

1.98 

1.05 

1.10 

0.86 

1.57 

0.76 

1.29 

1.0 

0.7 

2.09 

1.72 

0.93 

0.93 

1.16 

0.85 

0.86 

0.70 

1.31 

1.0 

0.6 

2.37 

1.28 

0.50 

0.78 

0.56 

0.79 

0.88 

0.66 

2.01 

1.0 

0.5 

2.73 

1.13 

0.23 

0.69 

0.31 

0.72 

1.14 

0.64 

1.41 

0.9 

0.9 

1.62 

2.15 

2.22 

1.07 

2.17 

0.77 

2.35 

0.66 

1.55 

0.9 

0.8 

1.84 

1.59 

2.08 

0.83 

1.92 

0.69 

1.55 

0.58 

0.85 

0.9 

0.7 

2.09 

1.44 

0.79 

0.81 

1.27 

0.64 

0.44 

0.55 

0.44 

0.9 

0.6 

2.37 

1.21 

0.86 

0.68 

0.58 

0.62 

1.29 

0.55 

0.70 

0.9 

0.5 

2.73 

1.00 

0.29 

0.59 

0.18 

0.63 

2.51 

0.50 

0.37 

0.8 

0.8 

1.80 

1.65 

1.62 

0.76 

1.37 

0.58 

1.55 

0.51 

1.91 

0.8 

0.7 

2.03 

1.59 

1.39 

0.76 

1.31 

0.56 

2.81 

0.45 

0.68 

0.8 

0.6 

2.30 

1.21 

0.56 

0.60 

0.40 

0.51 

0.62 

0.42 

0.62 

0.8 

0.5 

2.65 

0.96 

0.29 

0.52 

0.30 

0.49 

0.73 

0.41 

0.38 

0.7 

0.7 

2.00 

1.50 

0.84 

0.61 

0.80 

0.53 

1.49 

0.45 

1.07 

0.7 

0.6 

2.29 

1.10 

1.04 

0.53 

0.80 

0.45 

2.94 

0.34 

0.32 

0.7 

0.5 

2.61 

0.98 

0.33 

0.45 

0.11 

0.42 

0.32 

0.33 

0.32 

0.6 

0.6 

2.22 

0.93 

0.81 

0.40 

0.86 

0.48 

2.78 

0.31 

0.55 

0.6 

0.5 

2.52 

0.74 

0.48 

0.41 

0.36 

0.36 

0.32 

0.25 

0.32 

0.5 

0.5 

2.44 

0.85 

0.59 

0.37 

0.46 

0.35 

1.37 

0.23 

0.27 
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TABLE 2 

Summary of all calculated models 


Ml 

M 2 


Tmax 

p ( 7 max ) 

^ s,max 

p ( 7 s , max ) 

7 max,rem 

p ( 7 max , rem ) 

^ s , max,rem 

p (^ s , max , rem ) 

( Mq ) 

( Mq ) 

(10^ cm) 

( 10 ^ K ) 

(10*^ g cm-3) 

( 10 “’ K ) 

(10'’ g cm-3) 

( 10 '* K ) 

( lO *^ gcm-3) 

( lO ** K ) 

(10'^ gcm-3) 







Resolution 

= lOOA : M 0 -‘ 





1.1 

1.1 

1.25 


sn 

r 79 

33 S 


HZ 

1.25 

215 S 

1.1 

1.0 

1.45 

3.02 

6.71 

1.67 

3.05 

1.23 

2.67 

1.10 

3.11 

1.1 

0.9 

1.67 

2.54 

4.13 

1.46 

4.31 

1.17 

2.74 

1.01 

1.73 

1.1 

0.8 

1.86 

2.21 

1.92 

1.28 

1.67 

1.08 

2.12 

0.96 

2.96 

1.1 

0.7 

2.11 

1.82 

1.10 

1.09 

0.90 

0.98 

2.47 

0.90 

1.51 

1.1 

0.6 

2.39 

1.56 

0.42 

0.91 

0.65 

0.95 

1.82 

0.82 

1.24 

1.1 

0.5 

2.72 

1.42 

1.91 

0.77 

0.88 

0.92 

0.67 

0.81 

0.83 

1.0 

1.0 

1.43 

3.08 

4.80 

1.57 

4.17 

0.98 

3.05 

0.86 

2.57 

1.0 

0.9 

1.64 

2.10 

3.32 

1.17 

2.14 

0.93 

1.15 

0.80 

1.70 

1.0 

0.8 

1.87 

1.99 

2.49 

1.17 

1.07 

0.84 

0.73 

0.74 

0.97 

1.0 

0.7 

2.13 

1.64 

1.05 

0.98 

0.62 

0.81 

0.38 

0.68 

0.45 

1.0 

0.6 

2.40 

1.45 

1.23 

0.89 

0.76 

0.77 

0.98 

0.66 

1.16 

1.0 

0.5 

2.75 

1.30 

0.18 

0.70 

0.32 

0.75 

0.90 

0.62 

0.78 

0.9 

0.9 

1.64 

2.17 

3.21 

0.95 

2.34 

0.78 

2.92 

0.66 

2.09 

0.9 

0.8 

1.86 

1.96 

1.24 

1.07 

1.05 

0.71 

1.26 

0.60 

0.97 

0.9 

0.7 

2.08 

1.48 

1.51 

0.84 

0.56 

0.66 

0.40 

0.59 

0.38 

0.9 

0.6 

2.35 

1.27 

0.67 

0.73 

0.39 

0.65 

0.56 

0.58 

0.72 

0.9 

0.5 

2.68 

1.07 

0.31 

0.61 

0.18 

0.59 

0.48 

0.50 

0.50 

0.8 

0.8 

1.83 

1.65 

2.20 

0.83 

3.35 

0.59 

2.92 

0.54 

1.31 

0.8 

0.7 

2.07 

1.51 

0.89 

0.81 

0.82 

0.58 

0.79 

0.46 

0.97 

0.8 

0.6 

2.35 

1.27 

0.56 

0.63 

0.43 

0.57 

0.43 

0.45 

0.37 

0.8 

0.5 

2.69 

1.19 

0.37 

0.59 

0.19 

0.50 

2.74 

0.40 

0.51 

0.7 

0.7 

2.01 

1.44 

1.12 

0.80 

1.10 

0.49 

0.86 

0.42 

0.79 

0.7 

0.6 

2.27 

1.12 

0.45 

0.63 

0.39 

0.47 

0.60 

0.35 

0.46 

0.7 

0.5 

2.60 

0.84 

0.45 

0.49 

0.30 

0.46 

0.35 

0.34 

0.32 

0.6 

0.6 

2.25 

1.13 

0.76 

0.58 

0.67 

0.49 

2.61 

0.28 

0.43 

0.6 

0.5 

2.55 

0.81 

0.44 

0.44 

0.30 

0.32 

2.64 

0.25 

0.35 

0.5 

0.5 

2.47 

0.85 

0.52 

0.42 

0.47 

0.40 

1.89 

0.22 

0.32 






Resolution 

= 500 k M 0 -‘ 





1.1 

1.1 

1.25 

3.91 

5.84 

2.28 

5.27 

1.18 

3.10 

0.99 

2.90 

1.1 

1.0 

1.46 

3.50 

3.81 

2.52 

3.87 

1.18 

2.41 

1.02 

1.63 

1.1 

0.9 

1.67 

2.58 

4.24 

1.74 

3.96 

1.10 

1.48 

1.01 

1.11 

1.1 

0.8 

1.91 

2.19 

2.23 

1.43 

1.67 

1.04 

0.54 

0.94 

1.01 

1.1 

0.7 

2.16 

1.66 

1.50 

1.18 

0.69 

0.96 

0.99 

0.85 

0.91 

1.1 

0.6 

2.45 

1.68 

0.70 

1.00 

0.57 

0.95 

0.74 

0.82 

0.89 

1.1 

0.5 

2.81 

1.70 

0.76 

0.82 

0.22 

0.90 

0.63 

0.78 

0.61 

1.0 

1.0 

1.43 

3.50 

4.27 

1.61 

5.76 

0.92 

2.97 

0.80 

1.82 

1.0 

0.9 

1.64 

2.95 

3.78 

1.98 

1.81 

0.90 

1.89 

0.78 

1.77 

1.0 

0.8 

1.87 

2.26 

1.32 

1.54 

1.32 

0.87 

0.91 

0.77 

0.73 

1.0 

0.7 

2.11 

1.87 

0.93 

1.19 

1.13 

0.84 

1.06 

0.73 

0.33 

1.0 

0.6 

2.40 

1.31 

0.65 

0.69 

0.31 

0.79 

0.72 

0.66 

1.11 

1.0 

0.5 

2.74 

1.49 

0.58 

0.78 

0.15 

0.78 

0.74 

0.63 

0.51 

0.9 

0.9 

1.64 

2.77 

3.84 

1.36 

2.94 

0.79 

2.66 

0.56 

1.19 

0.9 

0.8 

1.86 

2.17 

1.26 

1.40 

0.77 

0.64 

0.87 

0.57 

1.03 

0.9 

0.7 

2.11 

1.72 

0.50 

1.05 

0.85 

0.64 

0.38 

0.57 

0.40 

0.9 

0.6 

2.40 

1.30 

0.43 

0.88 

0.36 

0.63 

0.51 

0.52 

0.46 

0.9 

0.5 

2.75 

1.28 

0.24 

0.70 

0.20 

0.59 

0.60 

0.50 

0.36 

0.8 

0.8 

1.82 

2.16 

2.08 

1.15 

1.71 

0.57 

1.86 

0.47 

0.44 

0.8 

0.7 

2.06 

1.74 

0.80 

0.96 

0.73 

0.50 

0.82 

0.43 

0.53 

0.8 

0.6 

2.34 

1.38 

0.75 

0.81 

0.46 

0.50 

0.54 

0.44 

0.24 

0.8 

0.5 

2.67 

1.04 

0.33 

0.54 

0.11 

0.50 

0.49 

0.41 

0.17 

0.7 

0.7 

2.05 

1.63 

1.38 

0.73 

0.93 

0.51 

1.71 

0.37 

0.59 

0.7 

0.6 

2.31 

1.11 

0.55 

0.57 

0.52 

0.48 

3.13 

0.33 

0.22 

0.7 

0.5 

2.65 

1.07 

0.25 

0.49 

0.19 

0.44 

2.57 

0.34 

0.15 

0.6 

0.6 

2.24 

1.45 

1.21 

0.79 

0.64 

0.50 

1.90 

0.27 

0.49 

0.6 

0.5 

2.55 

0.93 

0.52 

0.49 

0.16 

0.40 

1.77 

0.25 

0.10 

0.5 

0.5 

2.47 

1.09 

0.52 

0.56 

0.29 

0.41 

1.50 

0.21 

0.38 



